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ABSTRACT 

We investigate the dependence of QSO Lya absorption features on the temperature 
of the absorbing gas and on the amphtude of the underlying dark matter fluctuations. 
We use high-resolution hydrodynamic simulations in cold dark matter dominated cos- 
mological models. In models with a hotter intergalactic medium (IGM), the increased 
temperature enhances the pressure gradients between low and high density regions 
and this changes the spatial distribution and the velocity field of the gas. Combined 
with more thermal broadening, this leads to significantly wider absorption features in 
hotter models. Cosmological models with little small scale power also have broader ab- 
sorption features, because fluctuations on the scale of the Jeans length are still in the 
linear regime . Consequently, both the amplitude of dark matter fluctuations on small 
scales and thermal smoothing affect the flux decrement distribution in a similar way. 
However, the 6-parameter distribution of Voigt profile fits, obtained by deblending the 
absorption features into a sum of thermally broadened lines, is largely independent 
of the amount of small scale power, but does depend strongly on the IGM tempera- 
ture. The same is true for the two-point function of the fiux and for the flux power 
spectrum on small scales. These three flux statistics are thus sensitive probes of the 
temperature of the IGM. We compare them computed for our models and obtained 
from a HIRES spectrum of the quasar Q1422-I-231 and conclude that the IGM tem- 
perature at 2: ~ 3.25 is fairly high, Tq ^ 15000K. The flux decrement distribution of 
the observed spectrum is fitted well by that of a ACDM model with that temperature. 

Key words: cosmology: theory - intergalactic medium - hydrodynamics - large-scale 
structure of universe ~ quasars: absorption lines 



1 INTRODUCTION 

Neutral hydrogen in the intergalactic medium (IGM) pro- 
duces a 'forest' of Lya absorption lines bluewards of the 
Lya emission line in quasar spectra (Bahcall & Salpeter 
1965, Gunn & Peterson 1965). The general properties of 
these LyQ absorption lines are remarkably well reproduced 
in hydrodynamic simulations of cold dark matter (CDM) 
dominated cosmologies (Gen et al. 1994, Zhang, Anninos 
& Norman 1995, Miralda-Escude et al. 1996, Hernquist et 
al. 1996, Wadsley & Bond 1996, Zhang et al. 1997, Theuns et 
al. 1998). In these simulations, intergalactic gas settles into 
the potential wells of CDM dominated sheets and filaments, 
aptly termed the cosmic web (Bond, Kofman & Pogosyan 
1996). These large-scale density fluctuations in the inter- 
galactic gas are responsible for the weak absorption lines of 
column density Ahi < lO^'^'^cm"^ . 

In the simulations, the typical widths of the peaks in 
the neutral hydrogen density, which occur when a line-of- 
sight crosses a sheet or filament, depend on the amplitude of 
dark matter fluctuations on small scales and by the smooth- 



ing of these small scale fluctuations by gas pressure ('Jeans 
smoothing', e.g. Bryan et al. 1999). The widths of the ab- 
sorption troughs are then determined by differential Hubble 
flow accross the spatially extended absorber (Hernquist et 
al. 1996), peculair velocities (Hui 1999), and thermal broad- 
ening. 

Widths of absorption features have traditionally been 
characterized by the &-parameter distribution of Voigt pro- 
file fits, obtained by deblending the absorption features into 
a sum of discrete absorption lines. Recent high-resolution 
simulations have shown that these absorption lines in a stan- 
dard GDM cosmology have 6-parameters substantially nar- 
rower than observed in high-resolution Keck spectra (The- 
uns et al. 1998, Bryan et al. 1999). Since the 6-parameters 
have been demonstrated to be at least partially thermally 
broadened (Haehnelt & Steinmetz 1998, Theuns et al. 1998, 
Hui & Rutledge 1999, Schaye et al. 1999a), Theuns et al. sug- 
gested that an increased temperature of the IGM might re- 
solve the discrepancy. 

At the relevant redshifts (2 2: 5), the temperature 
of the low density IGM is determined by the interplay be- 
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tween adiabatic cooling and photoheating, which introduces 
a tight relation between density and temperature which can 
be well approximated by a power law (Hui & Gnedin 1997) . 
Schaye et al. (1999a) used numerical simulations to demon- 
strate that the observed cut-off in the Doppler parameter 
distribution and its dependence on column density is the re- 
sult of this p — T relation. The exact form of the relation 
depends on the history of reionization and on cosmological 
parameters (Miralda-Escude & Rees 1994, Hui & Gnedin 
1997). For a given reionization history, higher temperatures 
are obtained for models with a higher baryon density Q,bh^ 
and a smaller expansion rate of the Universe. Theuns et 
al. (1999) confirmed that in a model with high baryon den- 
sity D.bh'^ = 0.0250, low matter density fim = 0.3 and 
h = 0.65 the 6-distribution of the simulated absorption lines 
fits the observed distribution significantly better than for 
the colder standard CDM model. 

However, the value of Qbh^ required is higher than the 
current best estimates Qth^ = 0.019 (Buries & Tytler 1998) 
and it was suggested that additional heat sources neglected 
so far in the numerical simulations might contribute signifi- 
cantly to the heating rate. These include the photo electric 
heating by dust grains (Nath, Sethi & Shchekinov 1999) and 
Compton heating by the hard X-ray background (Madau 
& Efstathiou 1999). Alternatively, Abel & Haehnelt (1999) 
pointed out that the optically thin limit for the HeH pho- 
toheating rate usually assumed in these simulations might 
lead to an underestimate compared to the actual HeH pho- 
toheating rate by a factor of ~ 2 — 3. 

Here we use high resolution hydrodynamic simulations 
to investigate how the IGM temperature and the amplitude 
of dark matter fluctuations affect the distribution of the gas 
and influences the widths of absorption features. The ef- 
fect of different heating rates and the influence of the am- 
plitude of dark matter fluctuations on the gas distribution 
are discussed in Section 3. In Section 4 we investigate the 
relative importance of the three thermal smoothing effects, 
Jeans smoothing, peculiar velocities and thermal broaden- 
ing, and the amplitude of the dark matter fluctuations on 
the 6-parameter distribution obtained by Voigt proflle fit- 
ting. In Section 5 we discuss the effects on the one- and 
two-point function of the flux and on the flux power spec- 
trum. Section 6 contains our conclusions. As the preparation 
of this paper was in its final stages a preprint by Bryan & 
Machacek (1999) was circulated that investigates some of 
the same issues. 



2 SIMULATION PARAMETERS AND 
PREPARATION OF MOCK SPECTRA 

We have simulated several different models, characterized 
by their total matter density fim, the value of the cosmo- 
logical constant Oa, the rms of mass fiuctuations in spheres 
of 8h~^ Mpc, as and the heating rate. The parameters for 
these models are summarized in Table |^. All these mod- 
els have ilm + f^A = 1. For our reference model LI we 
take the currently favoured values h = 0.65 (see Freed- 
man et al. 1998), Slfc/i^ = 0.019 (Buries and Tytler 1998), 

* We write the present day Hubble constant as Hq = lOO/i km 
s-^Mpc-i 



Table 1. Models simulated 



Name 




h 


eg 




X-ray 


log To 


A 


L0.3 


0.3 


0.65 


0.9 


HM/3 




3.93 


1.07 


LI 


0.3 


0.65 


0.9 


HM 


no 


4.07 


0.92 


L2 


0.3 


0.65 


0.9 


2xHM 


no 


4.19 


0.78 


L3 


0.3 


0.65 


0.9 


3xHM 


no 


4.26 


0.68 


L10-.65 


0.3 


0.65 


0.65 


HM 


no 


4.07 


1.05 


LlcrA 


0.3 


0.65 


0.4 


HM 


no 


4.07 


1.25 


L3b 


0.3 


0.65 


0.9 


3xHM 


no 


4.26 


0.59 


LX 


0.3 


0.65 


0.9 


HM 


yes 


4.14 


0.81 


S3 


1.0 


0.50 


0.5 


3xHM 


no 


4.25 


0.78 



(ilm, ^a) ~ (0.3, 0.7) (Efstathiou et al. 1999 and references 
therein), erg — 0.9 (Eke, Cole & Frenk 1996) and a helium 
abundance of Y = 0.24 by mass. The IGM in this model is 
photoionized and photoheated by the UV-background from 
QSOs. The other models are variations on this set. Models 
L0.3, L2, L3, L3b and LX have the same normalization of 
the power spectrum as LI, but models Lla.65 and LloA 
have smaller erg. For the standard CDM model S3 we im- 
pose the normalization erg = 0.5 deduced from the abun- 
dance of galaxy clusters (Eke, Cole & Frenk 1996). We have 
used CMBFAST (Seljak & Zaldarriaga 1996) to compute 
the appropriate linear matter transfer functions. 

These simulations follow the evolution of a periodic, 
cubic region of the universe and are performed with a mod- 
ified version of HYDRA (Couchman et al. 1995). HYDRA 
combines hierarchical P3M gravity (Couchman 1991) with 
smoothed particle hydrodynamics (SPH, Lucy 1977, Gin- 
gold & Monaghan 1977). One of us (TT) has parallelized 
this code for the Origin 2000 at DAMTP, Cambridge. These 
simulations use 64^ particles of each species in a box of 
co-moving size 2.5/i~^Mpc, so the SPH particle masses are 
1.14 X lO'* (/i/O.65)"^M0 and the CDM particles are more 
massive by a factor 57cdm/J^6. The gravitational soften- 
ing lenght is 5 kpc. All our simulations were run with the 
same initial phases to minimize the effects of cosmic vari- 
ance when comparing the different models. We have also run 
an additional model (L3b) which simulates a larger region 
(5.0/i~^Mpc) with 128'^ particles of each type and thus the 
same numerical resolution as the other models. A simulation 
typically requires 1000-1500 steps to reach z = ?>. 

The detailed expressions for the heating and cooling 
rates as a function of temperature and ionizing fiux are taken 
from Gen (1992) with some minor modifications (Theuns et 
al. 1998). We take the photoheating rate in the optically thin 
limit from Haardt & Madau (1996). This heating rate, de- 
noted as 'HM' in the ene column of the Table, is imposed on 
models LI, L1(t.65 and Llcr.4. Models L2 and L3 are identi- 
cal to LI, except that we have multiplied the helium photo- 
heating rate by factors two and three respectively (keeping 
the ionization rate constant). Model L3b has the same heat- 
ing rate as model L3. In model L0.3 we have reduced the 
helium photoheating by a factor 3. Model LX is identical to 
model LI, except that we have included Compton heating 
by the X-ray background as a function of redshift z at a. rate 
(Madau & Efstathiou 1999) 

Hx(z) = 1.251 X 10"^' (1 + z)*-^ exp(-(z/5)^)ne , (1) 

where Ti.x is the heating rate in erg cm"^ s~^ and rie is 
number density of free electrons in cm~"^. Finally, model 
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S3 is a standard CDM model, with the appropriate Haardt 
& Madau (1996) helium heating rate increased by a factor 
3. 

For each model we compute spectra along 1200 random 
lines of sight through the simulation box at a given redshift 
(usually = 3) and then scale the ionizing background flux 
by the factor A in the analysis stage such as to give a cho- 
sen effective optical depth t^s ■ The factor A required to give 
Teff = 0.33 at 2 = 3 is given in Table 0. We process these 
simulated spectra using the following procedure, designed to 
introduce the same biases as present in observational data 
from the HIRES spectrograph on the Keck telescope. Each 
spectrum is flrst convolved with a Gaussian with full width 
at half maximum of FWHM = 8 km s~^and re-sampled onto 
pixels of width 3 km s~^. Photon and pixel noise is added 
such that the total signal-to-noise is 50. The 'continuum' of 
the spectra is then fitted with the method described in The- 
uns et al. (1998). 

We also compare fiux statistics for our simulated spectra 
to an observed HIRES spectrum of QSO 1422-1-231, kindly 
provided to us by W. Sargent and M. Ranch. The Lya for- 
est of this QSO extends over a significant redshift range, 
over which the effective optical depth evolves significantly. 
In addition, the signal-to-noise ratio varies as a function of 
wavelength and flux. We model these effects for this particu- 
lar spectrum using the following detailed procedure (see also 
Ranch et al. 1997). We begin by dividing the spectrum in two 
halves. For each half, we take the simulation output at the 
appropriate redshift and scale the background flux such as 
to match the observed mean effective optical depth. We con- 
volve the spectrum with a Gaussian with full width at half 
maximum of FWHM = 6.6 km s^^and re-sample it to pixels 
of the same size as the observed spectrum. We calculate the 
noise properties of the QSO spectrum as a function of flux 
and add Gaussian noise to the simulated spectra appropriate 
for the flux in each pixel. In the observed spectrum we ex- 
clude regions with identifled metal absorption lines and the 
region close to the emission redshift which might be influ- 
enced by the proximity effect. Below we will show that this 
more ellaborate procedure produces flux statistics in bet- 
ter agreement with the data than the simpler one described 
earlier. We use the same automated version of VPFIT (Car- 
swell et al. 1987) to fit Voigt profiles to both simulated and 
observed spectra (see Theuns et al. 1998 for more details). 



3 BROADENING OF QSO ABSORBERS 

In this section we illustrate the effect of different broad- 
ening mechanisms on the widths of density peaks along a 
given line-of-sight. In particular, we discuss the relative im- 
portance of the three different thermal effects, namely Jeans 
smoothing, pressure-induced peculiar velocities and thermal 
broadening and compare these to the effect of the amplitude 
of dark matter fluctuations. In the following sections we will 
investigate how the different broadening mechanisms affect 
a number of different flux statistics. 



3.1 Thermal smoothing 
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Figure 1. Temperature-density relation for model LI (points) 
and a power-law fit (least absolute deviation; full line). Corre- 
sponding fits for models L3 (short dashed) and model L0.3 (long 
dashed) are shown for comparison. The 'quantized' appearance 
of the data points close to the full line is due to the usage of an 
interpolation scheme to compute temperatures. 
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Figure 2. Volume weighted gas density probability distribution 
for models L0.3 (long dashed line), LI (full line) and L3 (short 
dashed line). Dark matter distributions are indistinguishable for 
these models and are shown as the dotted line. Voids in the hotter 
model are less empty because of increased Jeans smoothing, and 
so a smaller fraction of the volume is occupied by very low density 
gas. 
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Figure 3. Probability distribution for the divergence of the pe- 
culiar velocity for models L0.3 (long dashed line), LI (full line) 
and L3 (short dashed line). The larger pressure gradients in hot- 
ter models lead to more pronounced expansion of the gas out of 
potential wells. 



3.1.1 The spatial distribution and the velocity of the gas 

Figure ^ shows the temperature-density relations for the 
LCDM models where the Hell photo-heating is (0.3,1,3) 
times that in the optically thin limit, respectively. Over the 
interval shown, this relation is well approximated by a power 
law. The slope of this power law is determined by the reion- 
ization history (Hui & Gnedin 1997), which is the same for 
all models plotted here. The temperature To at the mean 
density for z = 3 is also given in Table As expected. To in- 
creases with the heating rate as To oc ei^^^j.^ (Miralda-Escude 
& Rees 1994). 

In our models, baryons represent less than 15 per cent 
of the mass and increasing their temperature does not affect 
the dark matter distribution very much. In contrast, the gas 
distribution itself changes significantly as the temperature 
is raised, due to the larger pressure gradients that arise be- 
tween high and low density regions. The increased pressure 
gradient forces the gas to settle into a configuration with 
a smaller density gradient in order to reach pressure equi- 
librium (VP/p — —S/cj)), leading to broader structures. This 
effect is illustrated in Figure ^ where we plot the probability 
distribution function (PDF) of the gas density for the three 
models. Hotter models have a smoother density distribution 
with less empty voids and less gas at high density (and hence 
a narrower PDF), while the PDF of the DM density is un- 
changed. (Note that the shift between the PDFs of gas and 
DM shows explicitly that our simulations resolve the Jeans 
length.) 

Another consequence of larger pressure gradients is that 
pressure-induced peculiar velocities in the gas become larger 
(see also Bryan et al. 1999). Figure ^ shows the PDF of the 
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Figure 4. Top panels show the density and peculiar velocity 
along a sight line picked from the simulation, for models L0.3 
(long dashed line), LI (full line) and L3 (short dashed line). The 
dark matter profiles for these models are identical and are shown 
dotted. Panel (c) shows the corresponding spectra. Panel (d) 
shows spectra for the case where peculiar velocities are not in- 
cluded. Finally, panel (e) shows the spectrum for model L0.3, but 
computed with the temperature-density relation of models L0.3 
(long dashed line), LI (full line) and L3 (short dashed line). In 
the hotter models, density peaks are broader, expansion velocities 
higher and the corresponding absorption features wider. 

divergence of the velocity field for the gas at an overdensity 
> 3. The distributions shift towards more positive values of 
div Up, i.e. a larger fraction of the gas at the center of sheets 
is expanding, for hotter models. Figure ^, to be discussed in 
more detail below, shows explicitly for a typical sight line 
how pressure-induced forces generate peculiar velocities that 
lead to a smoothing of the gas with respect to the underlying 
dark matter. 

The temperature also affects the normalization constant 
which determines the mean ffux level for a given ionizing 
background. If the gas distribution were to be independent 
of temperature, then this constant would simply scale as 
T""'^ due to the temperature dependence of the recombina- 
tion coefficient. There is, however, less gas locked in dense 
regions in the hotter models and consequently the tempera- 
ture dependence is slightly weaker. 
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3.1.2 The widths of QSO absorbers 

Figure ^ shows the density and velocity fields and the cor- 
responding Lya spectrum along a typical sight line for dif- 
ferent models. Panel (a) shows how a higher temperature 
leads to broadening of the density peaks in real space. The 
density peaks in the gas distribution are also significantly 
broader than in the underlying dark matter due to this 
Jeans smoothing. Panel (b) shows the corresponding pecu- 
liar velocity fields. The gas exhibits significant velocity bias: 
whereas the dark matter is collapsing towards density peaks, 
the gas is infalling well away from the central peak but is 
expanding near the density maxima (see also Bryan et al. 
1999). Photoheating is in the process of redistributing gas 
in the shallow potential wells where it collected at higher 
redshift when the temperature was lower. This effect is ob- 
viously stronger for models with a higher heating rate. 

The combination of thermal broadening, Jeans smooth- 
ing and pressure induced peculiar velocities makes absorp- 
tion features broader for hotter models, as shown in panel 
(c). The relative contribution of each of these effects is shown 
in the final two panels. Panel (d) shows the case where no 
peculiar velocities are included in the computation of the 
simulated spectra. Comparing panels (c) and (d) we see that 
both peculiar velocities and Jeans smoothing are important 
in determining the shape of the line wings. The wings of 
lines in model L3 are broader than those for model L0.3 
in panel (d) [without peculiar velocities). The difference is 
even greater in panel (c) (i.e. with peculiar velocities). Fi- 
nally, panel (e) shows spectra for a given density distribution 
and peculiar velocity field (model L0.3), but with different 
temperature-density relations imposed. These spectra which 
now only differ in the amount of thermal broadening look 
much more similar, but still the shape of the absorption pro- 
files near the peak are more rounded in the hotter models. 
Note that there are also some differences in the line wings. 

3.2 The effect of small scale power 

Hui & Rutledge (1999) derived a theoretical expression for 
the shape of the b-distribution that fits the observations 
fairly well. They made the Ansatz that QSO absorption lines 
arise from peaks in the density field. Their model explains 
the observed power law tail of large 6- values as being due to 
the frequent low-curvature fiuctuations whereas the sharp 
cut-off at low fo is a consequence of high curvature peaks 
being exponentially rare. The model also predicts that the 
widths of the absorption features should depend on the am- 
plitude of the underlying dark matter ffuctuations. Models 
with little small scale power are expected to have on average 
broader density peaks and hence broader absorption lines. 

Figures ^ and ^ demonstrate that the density peaks are 
indeed wider in the extreme low erg model L1(t.4 (og, = 0.4) 
than in model LI (us = 0.9, both models have almost 
identical temperatures). The collapse of the density fiuctu- 
ations around 100 and 200 km s~^ is significantly further 
advanced in model LI compared to model L1(t.4 and the col- 
lapsed structures are much less extended in velocity space. 
As a result the absorption features in model L1(t.4 are much 
broader than in LI. 

However, as is also shown in the figure, the number of 
Voigt components fitted to the resulting spectra by the Voigt 
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Figure 5. Same sight line as shown in figure ^ for the low fluctu- 
ation amplitude model Llcr.4. Top two panels show density and 
peculiar velocity, for the dark matter (dashed line) and gas (full 
line). The gas density for model LI is shown dotted in panel (a) 
for comparison. Panel (c) shows the corresponding fitted spec- 
trum (with noise added, wiggly line), the fit obtained using VP- 
FIT (full line), and the individual Voigt components (full lines 
with vertical line at maximum absorption), ofl^set vertically for 
clarity. The numbers in brackets denote the width of the fitted 
Voigt profile in km s"'^ for the lines with A^hi > 10^^ '* cm~^. 
Panel (d) shows the same as panel (c) for model LI. 

profile fitting routine VPFIT are different: the spectrum 
of L1(T.4 is fitted with five components while that of LI is 
fitted with only three components. It is therefore not obvious 
that standard Voigt profile fitting routines should require 
larger 6-parameters to fit the more extended density peaks 
in models with smaller amplitude of dark matter fiuctuations 
where the peaks often show significant substructure. We will 
show in detail in the next section that this is generally not 
the case. 



4 THE DOPPLER PARAMETER 
DISTRIBUTION 

We compute the Doppler parameters from Voigt profiles fits 
to absorption features. For both our simulated spectra and 
the observed spectrum of the Zcm = 3.62 QSO 1422-1-231, 
we use the same automated version of VPFIT . 
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Figure 6. Number of absorption lines per unit width db and 
per interval dv in Hubble velocity as a function of line width. 
Included are lines with column density in neutral hydrogen 
A^Hi > 10^^ cm~^ for which VPFIT gives an estimated error 
in b, berr < 0.25b. Panel (a) compares models with increased he- 
lium heating L1-L3. Panel (b) compares models L0.3 and L3 with 
(thick lines) and without (thin lines) peculiar velocities. Finally, 
panel (c) compares model L0.3 without peculiar velocities with 
the temperature-density relations of model L0.3, LI and L3 im- 
posed. Model L3, without peculiar velocities, is shown as the short 
dashed line. 

4.1 The global shift of the distribution with 
temperature 

The general shape of the 6-parameter distribution is rather 
skewed, with a slow fall of in the number of lines larger 
than the mode, but a rapid drop in the fraction of lines 
narrower than the mode. We will refer to the position at 
small b where the fraction of lines starts dropping rapidly as 
the 'cut-off' . (A more rigorous definition is possible (Schaye 
et al. 1999a) but here we will only compare fe-distributions 
between different models and so we do not need a more 
accurate definition of cut-off, since the differences amongst 
the models are relatively large anyway.) 

The effect of increased heating on the b-parameter dis- 



tribution is shown in figure ga. Increasing the helium heat- 
ing rate (models LI to L3) shifts the cut-off in narrow lines 
from < 20km s~^to > 25km s~^and decreases the fraction 
of lines at the peak. In addition to that, the hotter model L3 
has significantly more broad lines with b > 40km s^^than 
LI. Note that we have not normalized the distributions to 
demonstrate that the increase in broad lines is not due to a 
decrease in the number of narrow lines. As discussed in the 
previous section, the global shift is due to the combined ef- 
fect of Jeans smoothing, pressure- induced peculiar velocities 
and thermal broadening. 

The effect of peculiar velocities on the Doppler parame- 
ter distribution has been controversial and both a narrowing 
and broadening has been claimed. Figure ^ compares the 
line-widths distribution for models L0.3 and L3 with and 
without peculiar velocities. The effects of peculiar velocities 
are more pronounced for the hotter model L3. Peculiar ve- 
locities broaden the absorption lines and also significantly 
decrease the number of lines around the maximum in the 
distribution. This behaviour can be understood by looking 
again at figure ^, which shows that pressure-induced peculiar 
velocities tend to oppose the infall into sheets, and this effect 
is larger for the hotter model. Consequently, peculiar veloc- 
ities do not change the distribution of ^-parameters very 
much in the cold model while in the hot model the absorp- 
tion lines become broader. 

Note that these pressure induced peculiar velocities are 
not due to shocking of the gas: the higher temperature 
are purely a result of increased photo-heating. It is photo- 
heating that drives the higher density gas into its lower den- 
sity surroundings, which results in the changes in the density 
and velocity divergence PDFs shown earlier. The high den- 
sity gas collected at the bottom of the dark matter potential 
wells at higher redshifts, when the heating rates were lower 
and the gas was colder. 

To illustrate the effect of pure thermal broadening, we 
have computed spectra from the numerical simulations af- 
ter having artificially changed the temperature-density rela- 
tion of the gas. The resulting Doppler parameters for model 
L0.3, without peculiar velocity contributions, and with the 
temperature-density relations of model L0.3, LI and L3 im- 
posed are shown in Figure |^. Increasing the amount of ther- 
mal broadening shifts the cut-off at small b to larger values. 
The position of the cut-off is thus a measure of the tempera- 
ture of the gas. This is at the heart of the method developed 
by Schaye et al. (1999a), who used numerical simulations to 
demonstrate that the IGM temperature can be measured to 
high accuracy by comparing the lower cut-off of the Doppler 
parameter distribution of simulated spectra to that of the 
observed distribution. Comparing the full line with the dot- 
ted line shows that with increased thermal broadening the 
number of broad lines increases as well. 

The effect of Jeans smoothing can be seen by comparing 
the dotted with the short dashed lines. Both these models 
are with peculiar velocities set to zero and with the T ~ p 
relation of model L3, but the first uses the density distri- 
bution of L0.3 and the second the one of L3. The extra 
broadening of the fines is thus purely due to the widening 
of the absorbers in real space, due to the increased pressure 
gradients in the hotter model. The wider structures have an 
increased differential Hubble flow resulting in wider absorp- 
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tion lines. Note that the magnitude of this effect is smaller 
than that of thermal broadening. 



4.2 The effect of small scale power 



0.08 



0.06 - 



In section 3.2 we showed that models with less small scale 
power have significantly wider absorption features. The ef- 
fect of this widening on the &-distribution (i.e., the proba- 
bility distribution of lines as a function of their width, here 
for lines with Nm > 10^"^ cm~^) is shown in figure ^, which 
compares models LI and Llcr.4, which have the same T — p 
relation, but different values of ag. The two models have al- 
most identical 6-parameter distributions, which can be un- 
derstood as follows. When fitting a Voigt profile to a rela- 
tively broad absorption feature with a non- Voigt shape, as 
often occurs for models where the broadening is not ther- 
mally dominated, there is a trade-off between fitting a sin- 
gle broad Voigt line (which might not fit particularly well) 
or fitting multiple narrower lines. VPFIT was developed 
to 'deblend' absorption features into thermally broadened 
Voigt profiles and generally selects the option of fitting more 
but narrower lines with thermal shapes. Therefore, even for 
the wide absorption features in model LlaA the widths of 
the fitted Voigt profiles are still closely related to their "ther- 
mal" width. The 6-distribution for models LI and L1(t.4 are 
thus very similar, since they have the same temperature. 

We have fitted the two models with different as with 
another Voigt profile fitting routine, AUTOVP, kindly pro- 
vided to us by Romeel Dave (Dave et al. 1997) to check if 
this result is special to VPFIT . The 6-distribution obtained 
with AUTOVP is very similar to that obtained with VP- 
FIT , in particular, there is no strong dependence of the 
6-distribution on the amplitude of the dark matter fiuctu- 
ations. This result contrasts with the findings of Bryan & 
Machacek (1999), who did find a strong dependence of the 
6-distribution on erg using different profile fitting software. 

Finally, we confirm the result of Gnedin (1998) that 
changing the amplitude of dark matter fluctuations affects 
the column density distribution as shown in figure |^. At 
low column densities this is because many more weak lines 
are now required to fit the broader features. At the same 
time, LlaA has fewer strong lines and so the column density 
distribution steepens with decreasing ag. 



4.3 Comparison to the observed Doppler 
parameter distribution 

Figure |^ compares the Doppler parameter distribution of 
some of our models to the observed distribution for QSO 
1422+231. We have fitted the Lya -spectrum of this QSO 
using our automated version of VPFIT, thereby analysing 
observed and simulated spectra in the same way. We plot 
the 6-distribution of the fitted lines for the top and bottom 
half of the spectrum separately. 

Model LI reproduces the observed 6-distribution of the 
top half of QSO 1422-1-231 reasonably well, providing a good 
fit to both the cut-off at small 6, and to the tail of broader 
lines with 6 > 40 km s~^. However for the lower redshift 
half, LI seems to cut-off at too small values of 6, and in 
addition fails to reproduce the observed number of broader 
lines. Interestingly, the hotter model L2 does provide a good 



vpfit: 




20 40 
b (km/s) 

Figure 7. Fraction of lines with a given width, f(b), for models 
LI (full line) and Llo-.4 fitted with VPFIT (dotted fine) and 
using AUTOVP (dashed line). Although the structures in the 
extreme low erg model are significantly broader than in LI, this 
is not refiected in the 6-distribution. The b-distributions obtained 
using VPFIT and AUTOVP are very similar. 
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Figure 8. Column density distribution for model LI (open 
squares) compared to hlaA (filled squares). The low normaliza- 
tion model has significantly less high column density systems. 



fit to the bottom haff of QSO 1422+231 over the whole 6 
range. A more thorough investigation (Schaye et al. 1999b) 
demonstrates that the low redshift half of QSO 4122+231 is 
indeed hotter. 

Turning our attention to a comparison between the dif- 
ferent models, we see that models LI and LX produce very 
similar 6-distributions. Compton heating by the hard X-ray 
background is not able to raise the temperature of the gas 
probed by these lines significantly above the value which 
results from helium photoheating as computed by Haardt 
& Madau (1996). However there is a significant difference 
in temperature in regions of very low densities between LI 
and LX, which might be measurable with other flux statis- 
tics. Finally, models L2 and S3 also have very similar 6- 
distributions. Both models have a temperature of ^ 15000 K 
at mean density. In the standard CDM model S3 {ilm = 1), 
the adiabatic cooling time scale is smaller than in the ACDM 
model (r2„i — 0.3, SIa ~ 0.7). This explains why the Doppler 
parameter distribution are similar despite the larger Hell 
heating rate in model S3. 
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Figure 9. Fraction of lines with a given width, /(b), versus b 
for hnos with A'^jji > lO'^ cm~ with formal errors Oen- < 0.256, 
obtained from VPFIT . Symbols with (Poissonian) error bars are 
from fitting the spectrum of QSO 1422+231 using our automated 
version of VPFIT , which produced ~ 350 lines that met our 
selection criterion. Top (bottom) panel contains the top (bottom) 
half of the observed spectrum, the corresponding redshift ranges 
are indicated. The top panel compares the data to models LI 
and LX, the bottom panel to LI, L2 and S3. LI reproduces the 
observed distribution well for the top half of QSO 1422+231, but 
the bottom half seems to be fit better by L2. Models LI and LX 
are very similar, as are models L2 and S3. 

5 OTHER FLUX STATISTICS 

In this section we investigate the effect of the thermal broad- 
ening mechanism and that of the amplitude of dark mat- 
ter fluctuations on smaU scales on the the one- and two- 
point functions of the flux and on the flux power spectrum. 
These statistics are in principle easier to compute than the 
6-distribution but are likely to be affected by a number of ob- 
servational biases. For example, Rauch et al. (1997) showed 
that the shape of the one-point function (flux probability 
distribution) evolves strongly with redshift, as a consequence 
of the rapid evolution of the effective optical depth. There- 
fore, given that a QSO spectrum spans a relative large red- 
shift range, such evolution should be accounted for in a com- 
parison with simulations. We have followed the procedure 
described in detail in Section 2 to compare our simulations 
to the spectral statistics of QSO 1422+231. 

5.1 The one-point function of the flux 

The shape of the one-point function depends strongly on 
the assumed mean absorption, the continuum fit, and on 
the noise properties. For example, figure [l^a shows the dif- 



ferences between the top and bottom halves of the Lya spec- 
trum of QSO 1422+231 due to redshift evolution of the mean 
absorption, cosmic variance, and non-uniform signal to noise 
over the spectrum. 

The importance of using the correct noise properties 
is illustrated in figure |io|d where we compare the one-point 
function for a simulated spectrum with Gaussian noise inde- 
pendent of flux (dotted line) and with the more ellaborate 
method described in Section 2 (short dashed line). There 
are signiflcant differences, especially for saturated pixels 
{F = 0). Adding noise broadens the peak around nearly sat- 
urated pixels considerably compared to the case where the 
noise properties are similar to that in the observed spectra. 
In addition, a relatively large fraction of the pixels around 
F = comes from a few rare high column density systems. 
To test this we have removed the four highest column den- 
sity systems from the spectrum of Q1422 (column densities 
between 10^^'^ and 10^^'^ cm~^) and recomputed the one- 
point function. These two functions now differ significantly 
around F = 0. Since such strong systems are relatively rare, 
the shape of the one-point function round F — suffers from 
small number statistics. 

For the fuU curve in figure the spectra were calcu- 
lated as before but we have now also continuum fitted the 
simulated spectra as described in Theuns et al. (1998). Com- 
paring the short dashed line with the full line shows that, not 
surprisingly, continuum fltting has a big effect on the shape 
of the one-point function around F — 1. The smallness of 
our simulation box probably exacerbates this problem but 
it should be kept in mind that the continuum level of ob- 
served QSOs is often rather ill defined, especially at higher 
z. In the remainder of this section, we will present one-point 
functions computed using the noise properties of QSO 1422, 
and including continuum fitting. Finally, figure [lo|; compares 
models L3 and L3b: missing large scale power has a relatively 
significant effect on the one-point function, but this should 
not limit comparison between different models. 

Figure compares the one-point function for the cold 
model L0.3 to the hot model L3. The trend is as expected: 
the cold model has relatively more pixels round F ~ 0.1, 
since a larger fraction of the colder gas has collapsed to 
higher densities, but less absorption around 0.2 < F < 0.7, 
since the hotter model is denser in the regions surrounding 
dense peaks. This is due to the increased pressure forces in 
the hotter model driving more of the gas out of the centres 
of potential wells into the lower density regions. The middle 
panel shows that the effect of less small scale power is sim- 
ilar to that of higher temperatures: Ll(j.65 has fewer pixels 
round _F = but more absorption in the low density regions 
around 0.2 < F < 0.7. Because of the difference in power, 
gas in Lla.65 is less collapsed in high density regions, and 
so there is correspondingly more gas in lower density regions 
to boost the one-point function in the region 0.2 < F < 0.7. 
Machacek et al. (1999) also showed that the shape of the 
optical depth distribution is a sensitive probe of the am- 
plitude of density fluctuations on small scales. However, as 
we have just demonstrated, the effect of temperature on the 
flux decrement distribution and that of the fluctuation am- 
plitude of the DM density are largely degenerate. Finally, 
figure |ll| c shows that model L3b reproduces the observed 
one-point function reasonably well. 
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Figure 10. One-point function of the flux. Panel (a): comparison of top half of the Lya forest of QSO 1422+231 (dashed line) with 
the bottom half (full line). Cosmic variance, optical depth evolution and varying signal to noise are responsible for the differences. Panel 
(b); simulated spectra at the corresponding redshift, assuming Gaussian noise with effective signal-to-noise of 50 and no continuum 
fitting (dotted line) and using the noise properties of QSO 1422+231 without (short dashed line) and with (full line) continuum fitting. 
Continuum fitting affects the high end significantly, whereas the noise properties are important at low flux levels. Panel (c): comparison 
of models L3 (dashed line) and L3b (full line) illustrating the effect of missing large scale power. 




F 

Figure 11. Comparison between one-point functions. In panel (a) model L0.3 (full line) versus L3 (dashed line); in panel (b): model 
LI (full line) versus model Llo-.65 (dashed line); in panel (c) model L3b (dashed line) versus QSO 1422+231 (full line). Increasing 
temperature has the same effect as decreasing eg. Model L3b reproduces the observed one-point function reasonably well. 



5.2 The two-point function of the flux 

The two-point function of the flux, P{Fi, F2, Av), is the 
probability that two randomly chosen pixels separated by 
a velocity interval Av will have transmitted fluxes Fi and 
F2 (Miralda-Escude et al. 1997). A convenient way of visual- 
izing this function is by considering the mean flux difference 
for Fi in a given flux interval, that is, select a flux interval 
5Fi and consider the mean flux difference 



AF {v;SFi) 




(Fi - F2)P{Fi,F2,Av)dF2 } dFi/SFi (2) 



The shape of AF{v; SFi) contains information on the 
profile of the absorbers, or in the traditional language, on 
the widths of the lines. For large velocity differences, the cor- 
relation between the pixels is lost and the two-point function 
contains no extra information over the one-point function. 
In this limit, the mean flux difference will tend to the mean 
flux minus the mean flux of the pixels in the flux interval 
under consideration. Since A_F is an averaged quantity, the 
interpretation of the two-point function is not always obvi- 
ous. 

Figure |l2| compares the two-point function for sev- 
eral models and for QSO 1422 and for two flux intervals. 
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-0.05 < Fi < 0.05 (flux class 1) and 0.05 < Fi < 0.15 
(class 2). The temperature of the gas is clearly important: 
the two-point function shows that pixels are correlated over 
larger scales for hotter models, as expected given the broad- 
ening of absorption lines with increasing temperature. As 
with the one-point function, there is also a rather important 
effect of missing long waves in our small simulation boxes. 
The effect of small scale power on the other hand is sur- 
prisingly small, the two-point functions of models LI and 
Llcr.65 are almost indistinguishable. Only for model L1(t.4 
which has very much less small scale power are the differ- 
ences significant. This is mostly due the much lower number 
of high column density absorbers in this low normalization 
model. 

Model L3b reproduces the two-point function of QSO 
1422 up to 30km s~^, but then appears less correlated 
for larger velocity differences than the data. Unfortunately, 
just as for the one-point function, ffux class 1 is relatively 
strongly influenced by the presence of a few strong absorbers 
in the data (compare the dashed with the dotted line). The 
relatively good correspondence between model L3b and the 
data suggests that the temperature of the IGM at mean 
density is high, <; 15000K at this redshift. 

5.3 Flux power spectrum 

We begin by discussing the relation between the flux power 
spectrum and the underlying spectrum of density fluctua- 
tions. For a Gaussian random field, the power spectrum of 
density fluctuations along a sight line, Piu{k), is defined as 
the ensemble average of the square of the Fourier transform 
of the density p{v), 



Spik) 



p{v)exp{—ikv)dv / (p) , 



(3) 



where V is the length of the sight line, and p is the density 
as function of the Hubble velocity v. It is related to the three 
dimensional power spectrum of the field, P3D(|k|), through 
an integration (Kaiser & Peacock 1991) 



PsD{q)qdq . 



(4) 



This equation shows that PiD(fc) has contributions from all 
waves with q > k. Since Jeans smoothing introduces a small- 
scale cut-off in the power spectrum of the three dimensional 
density fleld, it influences Piu{k) for the gas on all scales. 

The observed flux is related to the gas density in redshift 
space by the relation 



F = exp[-^(pgas/ < Pgas >)"] 



(5) 



and the flux power spectrum Ppik) is defined as the ensem- 
ble average of the square of the Fourier transform of F{v), 



(6) 



Spik) — / F{v) exp{—ikv)dv . 



Peculiar velocities and thermal broadening contribute to 
the small scale cut-off of Ppik). The factor A (oc A) in 
equation ^ is a normalization constant that determines the 
mean flux, and acts as a 'bias' factor between density and 
flux. The exponent a depends on the slope of the p-T rela- 
tion. Ppik) differs from Pmik) due to the combined effects 



of Jeans smoothing, redshift space distortions and thermal 
broadening, and in addition due to the non-linear relation ^ 

Figure ^ compares k Pp {k) for several models and for 
QSO 1422-1-231. On small scales still above the resolution 
limit in the dark matter, the power spectrum dependence 
PsD ~ k~^ implies that kPio{k) depends only weakly on fe, 
highlighting that the observed sharp cut-off in k Ppik) is due 
to thermal smoothing. To compute these power spectra, we 
performed the following steps. The simulated spectra were 
calculated in the same way as we before, i.e. , scaling the 
models to the same mean flux, and adding noise with the 
same properties as that of the spectrum of QSO 1422. We 
have chosen a somewhat smaller pixel size of 2.2km s"Ho 
facilitate the use of a fast Fourier transform routine to calcu- 
late Pf (k) . To calculate the power spectrum of the observed 
spectrum, we re-sampled the spectrum to the same pixel size 
as the simulated spectra and divided it into stretches with 
lengths equal to that of our simulated spectra. We took care 
to use streches that have a flux close to one at the edges and 
in addition multiplied the flux with a windowing function, 
in order to be able to use 'periodic' boundary conditions in 
the FFT. We then computed the Fourier transform for each 
stretch separately and averaged the power per mode. To esti- 
mate the influence of missing large scale power in our small 
simulation boxes, we compare models L3 and L3b, which 
are identical expect for the size of the simulated region. As 
we show below, the differences in the flux power spectrum 
between these models is very small. 

On large scales where gas follows dark matter, equa- 
tions ^ and ^ can be combined to infer the three dimen- 
sional power spectrum from the flux power spectrum (Croft 
et al. 1997). In contrast, on small scales, logfc(km~^ s)> - 
1.4, thermal smoothing strongly suppresses the amplitude 
of Ppik). Panel (a) demonstrates how this cut-off moves 
to smaller wavenumbers for higher temperatures. The solid 
curve is for the cold model L0.3 while the dashed curve is 
for the hot model L3. The dotted curve is for the larger box 
L3b and there is a hint of increased power on the largest 
scales but the effect seems to be small. 

Panel (b) shows that there is little dependence of the 
small scale flux power spectrum on erg. This may seem sur- 
prising at first but is largely due to the fact that we have 
rescaled all simulated spectra to the observed mean flux 
(see Croft et al. 1997 for a discussion of this point). The 
dot-dashed line is for model LlaA, but imposing the same 
normalization A as for model LI, thereby increasing the ef- 
fective optical depth from ~ .33 to ~ .43. The power spec- 
trum for this model now differs significantly from the oth- 
ers. Consequently, scaling models to the same mean fiux re- 
moves most of the dependence of the fiux power spectrum on 
the fluctuation amplitude of the DM distribution on these 
small scales. Given that the shape of Ppik) depends strongly 
on temperature suggests it can be used as an independent 
method to measure the temperature of the IGM. 

Panel (c) compares L3b to the observed spectrum of 
QSO 1422. The fuU line shows kPpik) for L3b but without 
adding noise to the spectra. Comparison to the dashed line 
shows that the shape of Ppik) for logfc (km^^ s)>-0.5 is 
set by the noise properties. Imposing the noise of QSO 1422 
onto our hot model L3b gives a simulated power spectrum in 
excellent agreement with that of 1422. Therefore, both the 
two-point function and the flux power spectrum on small 
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Figure 12. Two-point functions for several models and for QSO 1422+231 for flux class 1 (—0.05 < Fi < 0.05) and flux class 2 
(0.05 < -Fi < 0.15). Panel (a): models L0.3 (full line), L3 (dashed line) and L3b (dotted line). Top curves are for class 1 and are off-set 
vertically by 0.3. Panel (b) is the same as (a) but for models Llo'.4 (full line), L1it.65 (dashed line) and LI (dotted line). Panel (c) is the 
same as panel (a) but for model L3b (full line), QSO 1422 with the four strongest lines removed (dashed line) and QSO 1422 (dotted 
line). Panel (a) shows that temperature and box size have an important effect on the two-point function, in contrast to the amount of 
small scale power (panel b). There are significant differences between the simulated and observed two-point functions, but these are at 
least partly due to missing large scale power in the simulations (panel a), and to strong absorbers in the observations (panel c). 
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Figure 13. Flux power spectrum for several models and for QSO 1422+231. Panel (a): models L0.3 (full line), L3 (dashed line) and L3b 
(dotted line); panel (b): models LI (full line). Lie. 65 (dashed line) and Llo-.4 (dotted line). Model Lie. 4 but with the same ionizing 
background as LI (same factor A, see table |l|) is shown with the dot-dashed line. Panel (c) shows models L3b without (full line) and with 
noise (dashed line), and QSO 1422 (dotted line). The bottom axes refer to wave number fe, the top axes to the corresponding Hubble 
velocity v = I-k jk. Panel (a) shows that temperature has an important influence on the shape of Pp{k), and that the effect of missing 
long waves is small. Panel (b) shows that ug has relatively little effect, if models are compared at the same mean flux. Panel (c) shows 
that the shape of Pp{k) for k > —0.7 is set by the noise, and that L3b fits the power spectrum of QSO 1422 on small scales very well. 



scales suggest that the IGM temperature at « ~ 3.25 is 
high, ^ 15000K. 



6 DISCUSSION AND CONCLUSIONS 

We have used numerical simulations to investigate what de- 
termines the widths of Lya absorbers in QSO spectra. These 



absorbers arise when a line of sight intersects a peak in the 
neutral hydrogen density, which itself is due to gas falling 
into the sheet-like and filamentary potential wells of the un- 
derlying dark matter distribution. The widths of the ab- 
sorbers are set by two very different processes. Firstly, the 
amount of small scale power determines the widths of the 
dark matter potential wells. Secondly, Jeans smoothing of 
the gas due to its high temperature ensures it is distributed 
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much smoother than the dark matter on smaU scales. Fi- 
nally, the corresponding absorption features are wider than 
the absorbers in real space due to peculiar velocities and 
thermal broadening. 

The strong dependence on the amount of small scale 
power is due to the fact that the linearly extrapolated fluc- 
tuation amplitude on the Jeans scale is of order unity for the 
models of interest. For this reason, the collapse of structures 
on scales comparable to the Jeans length is less advanced in 
models with a lower amplitude of DM fluctuations. 

The effect of temperature on the widths of peaks in the 
gas distribution can be understood as follows. Larger tem- 
perature gradients between high and low density regions, as 
occur in hotter models, lead to increased pressure gradients 
that push gas out into the wings of the underlying dark mat- 
ter potential wells. Consequently, less gas collapses into high 
density regions and voids in the density distribution become 
less underdense, making the peaks wider. This smoothing 
process is particularly efficient for photoionized gas because 
it obeys a temperature-density relation where the gas tem- 
perature increases with density, and this increases the pres- 
sure gradient over its isothermal value. 

We have further investigated what effect the different 
broadening mechanisms have on the resulting absorption 
spectra. We used the Doppler parameter distribution as ob- 
tained with the standard Voigt-profile fitting routines VP- 
FIT and AUTOVP, the one- and two-point functions of 
the flux, and the flux power spectrum to characterize the 
absorption spectra. 

The Doppler parameter distribution is mainly deter- 
mined by the combined effect of the three thermal smoothing 
mechanisms, thermal broadening, thermally driven peculiar 
velocities and Jeans smoothing. It shows little dependence 
on the amplitude of dark matter fluctuations. VPFIT and 
AUTOVP are therefore rather succesful in deblending ab- 
sorption features into thermally broadened Voigt proflles. 
Individual Voigt components often trace the substructure in 
broader density peaks. This is the reason why the Doppler 
parameter distribution shows little sensitivity to the over- 
all width of the density peaks, which themselves depend 
strongly on the amplitude of dark matter fluctuations. 

The flux two-point function and the flux power spec- 
trum on small scales are also mainly determined by ther- 
mal broadening and show likewise no strong dependence on 
the amplitude of dark matter fluctuations. This is at least 
partially due to the fact that we have chosen to scale our 
simulated spectra to match the mean observed flux. The 
one-point function of the flux, however, depends on temper- 
ature and normalization of the DM fluctuation amplitude in 
a similar way. This degeneracy arises because the PDF of 
the density depends mainly on the rms fluctuation ampli- 
tude of the gas, which itself depends on both the fluctuation 
amplitude of the DM density and on the scale on which the 
gas distribution is smoothed relative to the DM distribution. 

We conclude that the Doppler parameter distribution, 
the flux two-point correlation function and the flux power 
spectrum on small scales are sensitive probes of the temper- 
ature of the photoionized IGM, while the one-point function 
of the flux can be used to probe the fluctuation amplitude 
on small scales once the temperature is known. 

When comparing our simulated spectra with a high 
resolution spectrum of QSO 1422-1-231, obtained with the 



HIRES spectrograph on the Keck telescope, we flnd that 
the Doppler parameter distribution, the two-point correla- 
tion function and the flux power spectrum are best matched 
by models that have a temperature of the IGM of J; 15000K 
at mean density at z = 3.25. This temperature is ~ 50 per 
cent higher than in our reference model LI that used the 
heating rates calculated for the Haardt & Madau spectrum. 
(Note however, that in this paper we have only varied the 
amplitude of the p — T relation, and not its slope. This pre- 
vents us from making stronger statements about the value 
of To.) The flux one-point function depends both on tem- 
perature and the level of small scale fluctuations. A ACDM 
model with this temperature matches the observed one-point 
function reasonably well. Once an accurate measurement of 
To is available, it should be possible to determine erg accu- 
rately from the Lya forest. 
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